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Abstract 

The new symmetry axis of a well-behaved electromagnetic beam advanced in paper Physical 
Review A 78, 063831 (2008) is not purely a mathematical concept. The experimental result reported 
by Hosten and Kwiat in paper Science 319, 787 (2008) is shown to demonstrate the existence of 
this symmetry axis that is neither perpendicular nor parallel to the propagation axis. 
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Introduction. — Finite electromagnetic beams have now played very important roles in 
diverse areas of applications such as optical trapping and manipulation ll|, optical rotating 



optical guiding 



3|, optical data storage, and dark- field imaging But an incredible 



'act is that the theoretical description of a well-behaved finite beam has not been satisfactory 



n 



5I, y, 0, y, y, q, q, I2I, Q since the advent of masers and lasers Q, ll5|, [iGj. It is well 

known that the vectorial property of a plane wave is described by the polarization state. A 
useful concept is the Jones vector that consists of the two mutually orthogonal transverse 

n n 

components [17]. But for a finite beam, the polarization state is not a global property [7[. 
Instead it is local and changes on propagation. As a matter of fact, due to the longitudinal 
component [5|, ll3|j, the polarization state is not sufficient to describe the vectorial property 
of a finite beam. Recently, I advanced a new symmetry axis pL8|] represented by a unit 
vector I and found that it is a global property and, together with the Jones vector of the 
angular spectrum, can describe the vectorial property of a beam. The purpose of this 
paper is to show that this symmetry axis is not purely a mathematical concept. In fact, a 
perpendicular I to the propagation axis corresponds to the uniformly polarized beam in the 



paraxial approximation 



Id, 



6|, lZI], and a parallel I corresponds to the cylindrical vector beam 
13|. The conversion from uniformly polarized beams to cylindrical vector beams has 
been experimentally realized jiol. [lol. In this paper, I will show that a symmetry axis I that 
is neither perpendicular nor parallel to the propagation axis has been observed by Hosten 
and Kwiat 20|] in a recent experiment on the Imbert-Fedorov effect 2ll. |22|. 

Transverse displacement. — In order to see the connection of I with the Imbert-Fedorov ef- 
fect, let me first introduce the transverse effect the displacement of the beam's barycen- 
ter from the characteristic plane formed by I and the propagation axis, in the linear approx- 
imation. The electric vector F(x) of a monochromatic beam traveling in positive 2;-axis is 
expressed as an integral over the angular spectrum, 

1 



Ffx) 



271 



f (fca;, ky) exp(ik ■ 'x)dkxdk 



J/' 



(1) 



where kz = {k"^ — kl — k'^Y^'^, and the electric vector f of the angular spectrum is factorized 
as the following form, 

f(fcx, ky) = maf{kx, ky). (2) 
The mapping matrix m is so defined in terms of I and the wave vector k that it is normalized 



as m m = 1, where the superscript T stands for the transpose. As in Ref. 



18| , I is set to lie 



in the plane zox and to make an angle with the propagation axis, 1(0) = cos O+e^,. sin 0, 



where |0| < | is assumed. The normalized Jones vector 



a 



Olr. 



(3) 



is assumed to be the same to all the elements of the angular spectrum. The amplitude scalar 
f{kx, ky) of the angular spectrum is assumed here to include a phase factor, 



(4) 



in circular cylindrical system, where fo{kp) is square integrable and is sharply peaked at 
kp = 0, and / is an integer. The transverse displacement yb of the beam's barycenter is 
defined as the expectation of the y coordinate, 

/ jF^yFdxdy 



JjF^Fdxdy' 

where the superscript f stands for the conjugate transpose. Substituting Eqs. 
(jl]) into Eq. ([5]) and noticing that /o is an even function of ky, one get 



Vb 



J Ja'm 



(5) 
and 

(6) 



/ j\fo\^dk^dky 

It is noted that this displacement is independent of the phase factor exp in the scalar 
function (1 11). When |0| ^ 59, where 59 is half the divergence angle, m is linearly approxi- 



mated as 



18| 



^cot0 



m = sgn(0) 

Substituting Eqs. ([3]) and ([7]) into Eq. ([6]), we have 



^ cot ^ 



(7) 



a 



Vb 



cot 0, 



(8) 



where a = a^aa is the polarization ellipticity of the angular spectrum, and a 



-i 

1 



is the Pauli matrix. 

Imhert-Fedorov effect: change of I hy refraction. — Next I will show that the observed 



Imbert-Fedorov effect in Ref. 



20| demonstrates the change of I by the refraction. Consider 



the refraction at an interface between two different dielectric media, ni = 1 and n2 = 1.515, 
as is depicted in Fig. [1], where the laboratory reference frame is denoted by XYZ, the 
reference frame associated with the incident beam is denoted by xyz, the reference frame 
associated with the transmitted beam is denoted by x'y'z', 6q and 6'q are the incidence and 
refraction angles of the propagation axis, respectively. The symmetry axis I of the incident 
beam in Ref. 20| is perpendicular to the propagation axis. Let us start with an incident 



beam the perpendicular I of which lies in the incidence plane. 
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FIG. 1: (Color online) Reference frames xyz, x'y' z' , and XY Z associated, respectively, with the 
incident beam, the refracted beam, and the laboratory for the refraction at an interface between 
two different dielectric media, ni and n2- 



1. Description of the incident beam. In frame xyz, the wave vector of the angular 
spectrum in the linear approximation is given by k = ( A;^^ ky k )^, where the wave number 
in medium Ui is k = Uiko, ko = 2tt/\q, and Aq is the vacuum wavelength. With = — |, 
the linearly approximated mapping matrix becomes 



m 




1 



\ k k ) 

In frame XY Z , the wave vector is transformed into 

/ 



^XYZ 



A,f-^n)k 



k sin ^0 + fcx cos 
k cos 6*0 — kx sin 6*0 



(9) 



(10) 
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where 



^ cos i} — sin i} ^ 



Dy{d) = 10 

^ sin cos ■}} 

is the rotation matrix of the reference frame around the y axis by an angle i). The mapping 
matrix is correspondingly transformed into mxYZ = Dy{—do)^ and is given by 



^ cos ^0-% sin 00 sin 



mxYZ 







1 



V 



sin 6*0 - ^ cos 6^0 cos 6*0 



(11) 



2. Operations of the refraction to the incident angular spectrum. Refraction includes two 
kinds of operation to the incident angular spectrum. One is to alter the wave vector and 
the mapping matrix by rotation, and the other is to change the Jones vector. So the electric 
vector of the refracted angular spectrum in the frame XYZ can be assumed to be 



'-XYZ 



m 



XYZ 



a'f, 



(12) 



where a' is also normalized. Let us look at the former operation in frame XYZ. For an 
incident plane wave of wave vector ffTOl) . the incidence plane is formed by the wave vector 
(fTOj) and the Z axis. So the incidence angle 9 of this plane wave is determined by 

^XYZ -Gz „ K 



cos 9 



cos 9n 



\^XYz\ k 
Denoting by 9' the refraction angle, the Snell law gives 



sin 9o. 



■a' ■ a • /)/ , "-1^^ a 
sm.9 = — sm 9 = sm9Q -\ — cos o'o. 

n2 n2 k 

The unit vector n perpendicular to this incidence plane is given by 

ez X kxYz h, 



n 



\ez X kxYz\ ksm9o 
in the linear approximation. Denoting A^^ = 9' — 9, the direction of the refracted wave vector 
is obtained through rotating the incident wave vector around n by an angle A9 within the 
frame XYZ. The wave number is changed by the refraction into k' = ^k = n2/co- As a 
result, the refracted wave vector in the frame XYZ is given by 



XYZ = ^D-\A9)kxYz, 
III 



where 



D-\A9) 



cos A6 



k sin 60 



;i-cosA^) sinA^ ^ 



k sin do 

- sin AO 



k sin 6 



sinA6' 



k sin 00 



sin A6' 



cos AO 



J 



Correspondingly, the mapping matrix of the refracted plane wave is obtained through the 
same rotation and is given by 

"^'xYZ = D^^{A6)mxYZ- 



Now we look at the latter operation. Eq. ( jTT|) shows that the first and the second 
column vectors of the mapping matrix are p- and s-polarized, respectively, in the zeroth- 
order approximation. So the refraction converts the normalized Jones vector of the incident 



beam into Ta, where T 



tp 

u 



2 sin 9q cos 



sin(^o + ^0) cos(0o - ^0) 
is the Fresnel transmission coefficient of the p-polarization, and 



2 sin 6'^ cos 6q 



sin(^o + O'o) 

is the Fresnel transmission coefficient of the s-polarization. Introducing the normalized 



Jones vector, 



a' = ire . 




(13) 



where = (ItpCtpp + Its^sP)^''^ is the normalization coefficient, the amplitude scalar /' in 
the electric vector ([T2l) is given by 

/' = Nf. (14) 

Eq. ( fT3l) means that the polarization ellipticity of the refracted angular spectrum is different 
from that of the incident angular spectrum and turns out to be 



tptg 



(15) 



noticing that tp and tg are both real numbers. 
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3. Description of the refracted beam in frame x'y'z' . The wave vector of the refracted 
angular spectrum is transformed into 

= Dy{9Q)\i.XYZ 



h' 

\K'j 



( COS gp J„ \ 

COS dL ^ 



\ 



k' 



(16) 



The mapping matrix is transformed in the same way, m' = Dy{9Q)m'-xYZi ^^"^ expressed 
in terms of k' as follows, 



m 



1 

^yl COS Bq 

'T' ^ 



-COS 9L 



^yi cosSq^cosSq ^ 
T sin e'g 



sin 8'q 



\ 



x' 



y 



(17) 



Finally, we arrive at the electric vector of the refracted angular spectrum in the frame x'y'z' , 



f = m'a'f 



(18) 



where m', a', and /' are given by Eqs. ffTTI) . ffT^ . and ([HD, respectively. Comparison of Eq. 
( IT71) with Eq. (171) shows that the symmetry axis I' of the refracted beam lies in the plane 
z'ox' , that is to say, in the incidence plane. The angle B' between I' and the propagation 
axis is determined by 



cote 



, cos 9q — cos 9'q 
sin 6'n 



(19) 



and is no longer equal to — f . The dependence of 0' on the incidence angle 6^0 is shown in 




10 20 30 40 50 60 70 80 90 



(degrees) 



FIG. 2: Dependence of O' on the incidence angle Oq for ni = 1 and n2 = 1.515. 



Fig. [2] for ni = 1 and n2 = 1.515. 
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FIG. 3: Dependence of y'^ on the incidence angle for experimental parameters: ni = 1, n2 = 
1.515, and Aq = 633nm. 

According to Eq. ([8]), the displacement of the refracted beam's barycenter from the 
incidence plane is 

y^ = -^cote'. (20) 

This is nothing but the Imbert-Fedorov displacement. With ni = 1, n2 = 1.515, and 
Ao = 633nm, the dependence of on the incidence angle 6'o for a = ±1 is shown in Fig. [3l 
which is in quantitative agreement with the experimental data [2^. 

Indistinguishahleness of different perpendicular I — When Hosten and Kwiat performed 
their experiment, they did not realize the existence of I. So it was not ensured that the I 
of their incident beam lay in the incidence plane. In the following, I will show that incident 
beams of different perpendicular I are indistinguishable in the Imbert-Fedorov effect. 

For an arbitrary perpendicular symmetry axis, I = — e^; cos $ — e^^ sin $, the mapping 
matrix in the linear approximation is given by 

^ cos $ sin $ 

(21) 



m 



— sin $ 



sin $ 

— cos $ 



y ^ cos $ + sin $ cos $ - sin $ 



Inserting a unit 2x2 matrix / = Q Q into the right-hand side of Eq. ([2]), one has 



f = mQ^Q&f = m"a"f, 



(22) 



where 



Q 



cos $ — sin $ 
sin $ cos $ 



(23) 



8 



is a rotation matrix in the Jones-vector space, m" = mQ'^, and a" = Qa 




Eq. 



fl22|) means that the electric vector of the angular spectrum can be expressed either in terms 
of m together with Jones vector a or in terms of m" together with Jones vector a". With 
Eqs. (EH) and (1231), one has 

1 \ 

1 



m = — 



\ 



' k 



(24) 



k J 



which is the same as the mapping matrix 
the polarization ellipticity. 



At the same time, rotation Q does not change 



a" = a"^aa" = a^Q^aQa = a^aa = a. 



(25) 



It is thus clear that when the perpendicular I of the incident beam is rotated around its 
propagation axis, the linearly approximated mapping matrix and the polarization ellipticity 
can be regarded as remaining unchanged, resulting in the same Imbert-Fedorov displacement 

This work was supported in part by the National Natural Science Foundation of China 
(60877055 and 60806041), the Science and Technology Commission of Shanghai Munic- 
ipal (08JC14097 and 08QA14030), the Shanghai Educational Development Foundation 
(2007CG52), and the Shanghai Leading Academic Discipline Project (S30105). 



[1] A. Ashkin, IEEE J. Sek Top. Quantum Electron. 6, 841 (2000). 

[2] L. Paterson et aL, Science 292, 912 (2001). 

[3] X. Xu, K. Kim, W. Jhe, and N. Kwon, Phys. Rev. A 63, 063401 (2001). 

[4] D. P. Biss, K. S. Youngworth, and T. G. Brown, Appk Opt. 45, 470 (2006). 

[5] M. Lax, W. H. LouiseH, and W. B. McKnight, Phys. Rev. A 11, 1365 (1975). 

[6] L. W. Davis, Phys. Rev. A 19, 1177 (1979). 

[7] D. N. Pattanayak and G. P. Agrawak Phys. Rev. A 22, 1159 (1980). 

[8] L. W. Davis and G. Patsakos, Opt. Lett. 6, 22 (1981). 

[9] R. H. Jordan and D. G. Hall, Opt. Lett. 19, 427 (1994). 

[10] K. S. Youngworth and T. G. Brown, Opt. Express 7, 77 (2000). 

9 



[11] M. Onoda, S. Murakami, and N. Nagaosa, Phys. Rev. E 74, 066610 (2006). 
[12] K. Yu. Bliokh and Yu. P. Bliokh, Phys. Rev. E 75, 066609 (2007). 
[13] C.-F. Li, Opt. Lett. 32, 3543 (2007). 

[14] H. S. Green and E. Wolf, Proc. Phys. Soc. London Sec. A 66, 1129 (1953). 

[15] H. Kogelnik, Appl. Opt. 4, 1562 (1965). 

[16] H. Kogelnik and T. Li, Appl. Opt. 5, 1550 (1966). 

[17] R. C. Jones, J. Opt. Soc. Am. 31, 488 (1941). 

[18] C.-F. Li, Phys. Rev. A 78, 063831 (2008). 

[19] H. Ren, Y.-H. Lin, and S.-T. Wu, Appl. Phys. Lett. 89, 051114 (2006). 

[20] O. Hosten and P. Kwiat, Science 319, 787 (2008). 

[21] F. L Fedorov, Dokl. Akad. Nauk SSSR 105, 465 (1955). 

[22] C. Imbert, Phys. Rev. D 5, 787 (1972). 



10 



